arXiv:1509.00978v2 [hep-ph] 12Jul2016 


CNU-HEP-15-06 


Enhanced Charged Higgs Production through VF-Higgs Fusion 

in W-b Scattering 

Abdesslam Arhrib^’^, Kingman Cheung^’^’^, Jae Sik Lee^’^, and Chih-Ting Lu^ 

^ Departement de Mathematiques, Faculte des 
Sciences et Techniques, B.P 4 I 6 Tangier, Morocco 
^ Physics Division, National Center for Theoretical Sciences, Hsinchu, Taiwan 
^ Department of Physics, National Tsing Hua University, Hsinchu 300, Taiwan 
^ Division of Quantum Phases and Devices, School of Physics, 

Konkuk University, Seoul 143-701, Republic of Korea 
^ Department of Physics, Chonnam National University, 

300 Yonghong-dong, Buk-gu, Gwangju, 500-757, Republic of Korea 

(Dated: July 13, 2016) 

Abstract 

We study the associated production of a charged Higgs boson with a bottom quark and a light 
quark at the LHC via pp —)■ bj in the Two Higgs Doublet Models (2HDMs). Using the effective 
W approximation, we show that there is exact cancellation among various Feynman diagrams in 
high energy limit. This may imply that the production of charged Higgs can be significantly 
enhanced in the presence of large mass differences among the neutral Higgs bosons via -Higgs 
fusion in the pp —)■ bj process. Particularly, we emphasize the potential enhancement due to 

a light pseudoscalar boson A, which is still allowed by the current data by which we explicitly 
calculate the allowed regions in {Ma, tan/3) plane, and show that the production cross section can 
be as large as 0.1 pb for large tan/3. We also show that the transverse momentum distribution of 
the b quark can potentially distinguish the W^-A fusion diagram from the top diagram. Finally, 
we point out further enhancement when we go beyond the 2HDMs. 
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I. INTRODUCTION 


A new scalar boson h was discovered in the run I of LHC with 7 © 8 TeV energies in 2012 
DEI- The combined measurement of the mass of the boson performed by the ATLAS and 
CMS collaborations based on the data from h —)■ 77 and h —)■ ZZ —)■ 4/ channels is ruh = 
125.09 ± 0.21 (stat.) ± 0.11 (syst.) GeV [3]. Furthermore, the measured properties of the 
new particle are best described by the standard-model (SM) Higgs boson |11[5]. 

The mission of the new LHC run at 13 TeV (and later upgraded to 14 TeV) is two 
folds: the hrst task is the improvement of the scalar boson mass and scalar boson coupling 
measurements and the second one would be to hnd a clear hint of new physics. By performing 
accurate measurements of the scalar boson couplings to the SM particles would be helpful 
to determine if the Higgs-like particle is indeed the SM Higgs boson or a Higgs boson that 
belongs to a higher representation, such as models with extra Higgs doublets, extra triplets, 
or singlets. Most of higher Higgs representations with extra doublet or triplet Higgs helds 
predict in their spectrum one or more singly- or doubly-charged Higgs bosons. A discovery 
of such charged Higgs bosons would be an indisputable signal of new physics. 

In the two-Higgs-doublet models (2HDMs) or the minimal supersymmetric standard 
model (MSSM), the charged Higgs boson can be abundantly produced both at hadron and 
e'''e“ colliders. At hadron colliders, the charged Higgs boson can be produced through 
several channels: 

• Production from top decay. If the mass of the charged Higgs boson is smaller than 
mt — mb, the production of tt pairs provides an excellent source of the charged Higgs 
bosons. If kinematically allowed, one of the top and anti-top quarks, say the anti¬ 
top quark can decay into H~b, competing with the SM decay of t —?• W~b. This 
mechanism pp ^ tt ^ tbH~ can provide an important source of light charged Higgs 
bosons and offers a much cleaner signature than that of direct production. 

• Single charged Higgs production. The most important ones are gb — )■ tH~ and gg — )■ 

tbH~ [ 6 ]. These are QCD processes, and thus the cross sections are expected to be 
large. We can also have a single charged Higgs boson produced in association with 
a gauge boson via the loop process gg —)• or the tree level process bb —)■ 

m- Similarly, the single charged Higgs boson can be produced in association 
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with a Higgs boson: qq' —)■ W^* —)■ (j)H^ where 0 denotes one the the three neutral 
MSSM Higgs bosons [S]. Most of these processes are of the Drell-Yan type, they are 
expected to give substantial cross sections only for the charged Higgs mass below about 
200 GeV. 

• Single charged Higgs boson production associated with a bottom quark and a light 
quark qb —)■ q'H^b in the MSSM framework in which the neutral heavier Higgs bosons 
are almost degenerate j9]. 

• Charged Higgs pair production through qq annihilation HDl or gluon fusion. 

• Resonant charged Higgs production cs cb —)■ [TT] . 

At the Tevatron and LHC, detection of light charged Higgs boson with Mjj± < rrit — mb 
is straightforward from ti production followed by the decay i —?■ bH~ or t bH~^. Such a 
light charged Higgs boson can be detected for any value of tan/? in the rz/ decay which is 
indeed the dominant decay mode. The ATLAS and CMS have already had an exclusion on 
bH~^) X B{H^ —)■ ru) based on this decay channel [121 113 ] • 

In the MSSM and 2HDMs, the heavy charged Higgs boson with Mh± ^ mt would decay 
predominantly into tb. The experimental search is rather difficult due to large irreducible 
and reducible backgrounds associated with —)■ tb decay. However, in Refs. HI it has 

been demonstrated that the —>■ tb signature can lead to a visible signal at the LHC 

provided that the charged Higgs mass is below 600 GeV and tan/3 is either below < 1.5 or 
above > 40. An alternative decay mode to detect a heavy charged Higgs boson is —>■ rz/ 

[T^ . even if such a decay is suppressed for heavy charged Higgs bosons, it has the advantage 
of being much cleaner than —)■ tb. Recently, a new technique using the jet substructure 

for the heavy charged Higgs boson decaying to tb has been proposed in |16j . 

In the MSSM, the branching ratio of the decay mode B{H^ — W^h) could at best be at 
the level of 10% for low tan/3 while in the 2HDM-l|^it could dominate over B{H^ —)■ tb). 
Therefore, —)■ W^h could be an alternative channel to discover the heavy charged Higgs 

boson at the LHC HT]. Similarly, when the CP-odd Higgs boson A is light enough, the decay 
of W^A could be the dominant one in the 2HDMs and could also be used to search 

for heavy charged Higgs bosons. Finally, in the models with higher Higgs representations 

* See Section irri for classification of 2HDMs. 
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such as the triplet representation of the Higgs, the charged Higgs boson could decay into 
W^Z with a signihcant branching fraction [TH]. This decay channel could lead to isolated 
leptons in the hnal state and could be used to distinguish between models with charged 
Higgs bosons. 

The aim of this work is to study singly-charged Higgs boson production in association 
with a bottom quark and a jet g' with the subprocess qb — q'H^b. Such a process had 
been studied for the hrst time in Ref. |9] which showed that the rate is rather small in the 
MSSM due to a huge cancellation between the top- and Higgs-mediated diagrams as we will 
show. In the present study, we discuss the production rate of this process and its sensitivity 
to tan in the 2HDMs where the masses of the heavier Higgs bosons are not hxed by one 
mass parameter as in the MSSM. Specihcally, we demonstrate that the process possesses 
destructive interference between the s- and f-channel diagrams, which signihcantly reduces 
the cross section. Especially, when the two heavier neutral Higgs bosons are decoupled from 
the lightest one and they are degenerate, the cross section is canceled to a large extent. In 
addition, we show that with a relatively light CP-odd Higgs boson, which is still allowed 
by the current data, the production cross section of the charged Higgs boson via hP^-Higgs 
fusion in the pp —)■ bj process can be signihcantly enhanced at the LHC. 

The organization of the work is as follows. In the next section, we write down the 
framework for the 2HDMs, provide analytic understanding of the process in terms of the 
2 —)■ 2 subprocess, and also describe the full 2 —)■ 3 process in detail. We present the 


numerical results in Sec. HI Some cases beyond the 2HDMs are considered in Sec. IV and 
we conclude in Sec. El 


II. qb q'H+b IN TWO HIGGS DOUBLET MODELS 
A. Brief review of two-Higgs-doublet models 

In 2HDMs the electroweak symmetry breaking is performed by two scalar helds <hi and 
$2 which are parameterized by : 


$1 = 


^ (ui + 0? + iai) 


$2 = 


{V 2 + (j )2 + 10-2) 


( 1 ) 


t For an overview, see Ref. m- 
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TABLE I. Classification of 2HDMs satisfying the Glashow-Weinberg condition [22] which guaran¬ 
tees the absence of tree-level FCNC. 



2HDM I 

2HDM II 

2HDM III 

2HDM IV 

vf 

0 

1 

0 

1 

di 

1 

0 

1 

0 

hi 

0 

1 

1 

0 

hi 

1 

0 

0 

1 


We denote Vi = vcos/3 = vcp and ^2 = nsin/3 = vsp. The parameterization of the general 
scalar potential which is gauge invariant and possesses a general CP structure can be found 
in PD]. In the present study we are mainly interested in Higgs coupling to fermions and 
gauge couplings to be listed slightly later. 

The general structure for Yukawa couplings is given in the following interactions 

- Cy = Ku^Q'^ (fra) $2 + hd (Ir Q'^ (fra) $ 2 ) 

+ hiluL^ {iT 2 ) $ 2 ) + h.c. ( 2 ) 

where = {ur , d^), L'^ = {ur , Ir), and d)* = frad** with 



We note that there is a freedom to redehne the two linear combinations of $2 and $1 to 
eliminate the coupling of the up-type quarks to <hi PH. The 2HDMs are classihed according 
to the values of 2 Vi ,2 Table 

To define the Higgs mass eigenstates, we first rotate the imaginary components a* and 
the charged ones and 0 ^ in order to obtain the would-be-goldstones and that 
would be eaten by the longitudinal components of the Z and bosons. These rotations 
result in an CP-odd state a = A = —S/jOi -|- 0/302 and a pair of charged Higgs bosons 
= —S/ 30 ^ -b Cj 3 (j)f. In the most general case with CP violation, the mass eigenstates 
of the neutral Higgs bosons are obtained by diagonalizing the 3x3 mass matrix ATq by 
an orthogonal 3x3 mixing matrix O that relates the interaction eigenstates to the mass 
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eigenstates as follow: 


{cl)l4>2,a)a=OUHi,H2,Hs) 


(4) 


such that O^M^O = dia.g{Mfj_^, Mfj^, M]j^) with the ordering of Mh^ < Mh 2 < Mh^- Here 
the states Hi do not have to carry any dehnite CP-parity and they have both CP-even and 
CP-odd components. 

After identifying the Yukawa couplings by 

_ V2mu 1 , _ V^rud 1 , _ V^mi 1 

r^u -; n^d- - d -^— d — ’ ~ - 1 -^— i — ’ 

V Sf} V r]lcp + rjpfi v 7]{cf} + 

one can easily obtain, from the above Lagrangian, the following Higgs-fermion-fermion in¬ 
teractions 


m ,. 


= — 


V 


V 


u 


O 


4>2i 


S/3 


^ —Oai 75 
5/3 


u 


Hi 


T (vfo^ii + ■ hfS/3 - 72^/3 ^ 

d I j j ^ d \ d ^0,1 75 

hi C/3 + h2S/3 


V hlC/3 + h2S/3 


H 


J (+ V20^2i • his/3 - h2C/3 ^ 

1 III I I I 

L V hlC/3 + h2S/3 hlC/3 + h2S/3 


H 


( 6 ) 


and 


^H^ud 


V^m / r]{sf} - T]iCj3 \ 
V Vhic/3 + h2S/3/ 


(7) 


where Pl,r = (1 Y 7 ^)/ 2 . 

Before moving to the next subsection, we present the mixing matrix O in the CP- 
conserving case in terms of the mixing angle a. In our numerical study, to deliver our 
hndings more clearly, we focus on the CP-conserving case. In this case the matrix O takes 
the following form: 


O 


( 


— sin a cos a 0 
cos a sin a 0 


( 8 ) 


V » ° V 

assuming i /3 is the pure CP-odd state or i /3 = A. In this notation, the decoupling limit of 
the 2HDM |23], which seems to be favored by the current LHC data, is /i — a —)■ 7 r/ 2 : 
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FIG. 1. Feynman diagrams for 2 —)• 2 subprocesses: W~^b — H'^b (a) and (b), W~^b — H'^b (c) 
and (d). 


Ofp^i = — sin a —)■ cos j3 , = cos a —)■ sin (5 ; 

0(j,2i = cos a —)■ sin (3 , 0,^22 = sin a —)■ — cos /? . (9) 


B. Subprocess W~^ b —> b and unitarity 

In this subsection, we present the amplitude of the process qb ^ q'H^ b in the effective W 
approximation. In this process, the dominant contribution comes from the region where the 
W boson emitted from the incoming quark q is close to on shell and one can approximately 
represent the process by the W boson scattering with the incoming b quark or anti-6 quark 
to give H^b or H^b in the hnal state: 

W+{qi) b{pi) ^ H+{q2) b{p2). 

W+{qi) b{pi) ^ H+{q 2 ) b{p 2 ). (10) 

The process W^b — H~^b receives contributions from Fig. [^a) a f-channel diagram with 
the neutral Hi exchanges and Fig. [^6) a s-channel diagram with top exchange. While the 
process W~^b —)■ H^b receives contributions from Fig. [^c) a f-channel diagram with the 
neutral Hi exchanges and Fig. [^d) a w-channel diagram with top exchange. The relevant 
interactions needed for these two subprocesses can be obtained from the Yukawa interactions 
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given by Eqs. 0 and Q and from the covariant derivatives: 

^Hibb = 


where 


^H±tb — d 


b {gf + ig[ 75 ) bHi, 

b {cl Pl + cr Pr) tH~ + h.c., 


2mw 

gmb j 


y/2mw 

Pw±tb = -gjV2 {i'f^.PL b) + h.c., 


PuiH^w^ — (>S'i + iPi) 


H- it d^] Hi 


hh+^ + h.c., 


Cf}0^2i , Pi Oai , 


(11) 


( 12 ) 


and 


n ’rrit 1 

Cl = tan , cr = — 


O, 


rrib tan/? ’ 


gi = 


cp 


, gf = -iwil30a 


in types II and IV and 


Cl = 


tan/? ’ 


Cr = 


rrit 


rrib tan/? ’ 


gf = 


O, 


4>2i 


SfS 


p Odi 

gi = 


tan /? 


in types I and III. 

The amplitude of each diagram for W^{qi) h{pi) —)■ id+(g2) ^(^2) reads 


= __ 

Amw{t — Mjj) 


( 13 ) 


( 14 ) 


{Si + iPi) {q2 + PHTc^{qi) [u{p2) (^f + igf 75 ) u{pi)\ , 


■^ib) = -o - 2\ [(^Lu{p2)llt(f{qi)PLu{pi) + CRmtu{p2)(f{qi)PLu{pi)] , ( 15 ) 

lrrt\Y\s — rrti) 

where s = {pi + qif = (p2 + ?2)^ t = {P1-P2Y = (?2 -gl)^ and u = {pi-q2f = {P2-qif 
and e^{qi) denotes the polarization vector of W~^ boson. The amplitudes for the (c) and 
(d) diagrams in Fig. [^can be obtained by replacing n(pi,2) with ^(pig) and (s — mf) with 
(n-mf). 

In the high-energy limit, s, |t|, |m| 3 > , Mjj,, we hnd that 


M 


(a)+(b) 


'£m^)+m 


^ s'""'-'’ 

~ 4m^ 


2Ci + 5^(S.9f + P.9f 


- Pp) + + Pg?) 

i 

+ PPigf) u{p2)Plu{pi, 


u{p2)Pru{pi) 
( 16 ) 


where we have taken the longitudinally polarized W or e^(gi) ~ qi/irtw'- <?i = Pt ~ Pi with 
Pt = s for the diagram ( 6 ) and qi = Pt + P2 with pj = u for the diagram (d), respectively. 



















denoting the four-momenta of the exchanging top quark with pt. Incidentally, the square 
of the 4-momenta of the internal neutral Higgs is = {pi —P 2 Y — We note that the cr 
term, which is suppressed by rrit/y/s, is neglected here. In types II and IV, the cr term could 
be important when tan f3 < ~ 7. As shall be seen, the total cross section takes its 

smallest value at tan ~ 7. When tan /9 > 7, compared to the cr term, the cr term could 
be safely neglected when y/s/rrit 3> {mt/rrib)/ tan^/3. On the other hand, in types I and 
III, the Cr term can be neglected only if \/s/mt S> mt/rrib. Therefore, the high-enegy limit 
should be applied with more cautions at the LHC for types I and HI. But, for the 2HDM 
types I and HI, the production cross sections are suppressed by 1/ t&v? (3 with increasing 
tan (3 and the largest value with tan /3 = 1 is only ~ 30 fb, as shall be shown. 

The amplitude M(a)+(b) for the 6-initiated processes consists of the contributions from 
the f-channel Higgs-exchange diagrams (a) and the s-channel top-exchange diagram (b) . 
The Cr term in the second line is from the s-channel diagram and all the others from the 
t-channel ones. Therefore, the high-energy limit has been obtained by taking s/(s — ml) k. 
tj {t — MffJ 1. On the other hand, the high-energy limit of the amplitude A4.[c)+{d) for the 
^initiated processes can be obtained by replacing m(pi, 2 ) with f (pig) in Eq. ( [l^ and taking 
u/{u — mD ^ tj— ^ 1. We note that the high-energy behavior of Ai(^a)+ib) is different 

from that of M.(^c)+(d) especially when s is not large enough and there are non-negligible parts 
of phase space in which the high-energy limits u/{u — ml) ~ 1 in A4(c)+(d) are 

much more difficult to achieve than the corresponding ones s/(s — ml) k, t/{t — Mjj) k. 1 


in Otherwise, the expression given by Eq. (16) can be applicable for both the b- 

and 6-initiated processes. 


The high-energy limit expression Eq. (16) contains two non-interfering terms both of 
which grow as a/^ and therefore the absence of these unitarity-breaking terms require the 
following three types of sum rules: 


2cR + Y.^S,gf + P,gf) = Q, 

i 

E S’sf = E p-sf , 

i i 

E s.sf = E ^<9? = 0 ■ ( 1 ^) 

i i 

The hrst one gives the relation between the charged Higgs coupling to t and b quarks {cr) 
and the sum over the Higgs states of the scalar and pseudoscalar products {gfSi + gf Pi) of 
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the neutral Higgs couplings to b quarks and those to the charged Higgs and W. The second 
relation shows the sum over the Higgs states of the scalar products should be the same as 
that of the pseudoscalar ones. And the third relation implies that there is no CP violation 
if the scalar-pseudoscalar products are summed over the three Higgs states. 

These interesting sum rules can be explicitly checked in each 2HDM. In types H and IV, 
using the orthogonality of the mixing matrix O, we hnd that 

- tan/^Oj^J = - tan/3 , 

i i 

= - tan/3 ^ = - tan/3 , 

i i 

^ S.jf = 5/P.jf = 0 . (18) 

i i 

With Cl = tan /3, the unitarity conditions are satished automatically. On the other hand, in 
types I and HI, we hnd that 

= l/tan/3 , 

i i 

i i 

= Vasf = •)■ (19) 

i i 

With Cl = —1/ tan/3, the unitarity conditions are again satished automatically. 

This is the proof for the unitarity of the subprocess W~^b — bH~^ in the high energy limit 
in the general 2HDMs with or without CP violation. The same proof also applies to the 
case of b initiated subprocess W~^b — bH^. 

C. The full process qb —)■ q’H^b 

After discussing the essence of the physics involved in the 2 —)■ 2 subprocess, we shall 
describe the full 2 —3 process 0 We shall consider the CP-conserving case for simplicity, 
unless stated otherwise. In this case, without loss of generality, we identify Hi = h, H 2 = 
H, and = A, where h and H denote the lighter and heavier CP-even Higgs bosons, 

t For a full consideration of NLO corrections, one may need to take account of the 2 —>■ 4 process: 
qg —>■ qH~^bb. We leave this part for further work. 
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FIG. 2. Feynman diagrams for qb —)• q'H'^b (a) and (b), qb — q'H'^b (c) and (d) where {q-,q') = 
(m, d), (c, s) and Hi = h, H, A. The processes with {q, q') = (d, rZ), (s, c) are understood. 


respectively, and A the CP-odd one. The Feynman diagrams for the subprocesses qb -A 
q'H^b and qb -A q'H^b are shown in Fig. We stress at this level one important difference 
between the bottom-initiated diagram in Fig.j^b) and anti-bottom-initiated one in Fig.j^d) 
is that the former has a s-channel exchange top propagator while the latter has a u-channel 
one. Similarly, the fermion-line direction of the q can be reversed to include q ^ q' transition. 
Therefore, we have a number of initial states for production of (u, c, d, s) ® (6, b). We 
can then take the charge conjugate to obtain the H~ processes. 

The diagram in Fig. [^b) represents a top-induced process. If Mh± < rrit — mb, the top 
quark is produced on-shell, then followed by its decay into bH~^. This diagram is entirely 
dominant over the other diagrams. However, when Mh± > mt — the top quark is off- 
shell, and thus other diagrams also make signihcant contributions. In other diagrams, the 
charged Higgs boson appears being produced by WHi fusion, where Hi = h,H,A in the 
CP-conserving case. The coupling in the vertex W^H~Hi is a gauge coupling proportional 
to g and some mixing angles of the Higgs sector, and independent of different types of 
2HDMs. On the other hand, the dependence on the type of 2HDMs comes from the Yukawa 
couplings of Hi to b quark and the charged Higgs boson to tb. We list the relevant Yukawa 
couplings for 2HDMs from type I to IV in Table up to some normalizations. Incidentally, 
the non-vanishing neutral Higgs couplings to charged Higgs and W are given by 

Sb cos^/3 crj, 

82 = Sh = C/ 9 O 022 - S/ 3 O 012 = - sin(/3 - a ), 

Fs = Pa = Oa3 = 1, (20) 

using the form of O given by Eq. Q. 
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TABLE II. The bottom quark Yukawa couplings for h, H, A and that of the charged Higgs boson 
for 2HDMs of type I, II, III, and IV. The common factor for neutral Higgs boson is grrib/V^Mw 
while that for the charged Higgs is g/y/2M\Y- The chiral projection operators are Pl,r = (1T7'^)/2. 



Type I, HI 

Type H, IV 

hbb 

COS a 

sin 13 

sin a 
cos (3 

Hbb 

sin a 
sin j3 

cos a 
cos (3 

Abb 

+ cot (3 

— tan j3 

H-tb 

mb p 1 mt p 
tan /3 ^ tan /3 ^ 

mb tan 13 Pl + ^^^^Pr 


In the decoupling limit, we have cos{/3 — a) = 0 and sin(/3 —a) = 1. The contribution from 
the light Higgs h diagram is automatically zero because Sh = 0. The contributions from H 
and A are the same up to the 7 ® factor in the vertex if they are degenerate. If we look 
at the diagram more closely, the whole process can be regarded as Wh and Wh annihilation, 
as 2 —)■ 2 processes. It is easy to see from Fig. that for the Wh — H^b subprocess we 
have three t-channel diagrams with Hi = h,H, A in Fig. [^a) and one s-channel diagram 
mediated by the top quark in Fig. [^b). Similarly, for Wh —>■ H^h subprocess we have three 
f-channel diagrams with Hi = h,H,A in Fig. [^c) and one w-channel diagram mediated by 
the top quark in Fig. [^d). We have shown in the previous subsection using the effective 
W approximation that there is strong cancellation among the diagrams, and indeed all four 
diagrams will exactly cancel one another in the high energy limit. Therefore, if we employ 
a much lighter CP-odd Higgs boson, which is still allowed by the current data, we expect 
a strong enhancement to the production cross section of this process. Experimentally, one 
can use this process to search for the charged Higgs boson and investigate the effects of light 
CP-odd Higgs boson. Perhaps, a negative search would close out the entire window of light 
CP-odd Higgs boson. 


III. NUMERICAL RESULTS 

In this section, we hrst present some numerical results for the subprocesses W~^b —)■ bH~^ 
and W~^b —)■ bH^ for a given value of center-of-mass energy \/S and then consider the full 
process pp —)■ H'^bj in the 2HDM of type I (HI) and H (IV). 
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A. The 2 —)• 2 subprocess in the effective W approximation 

We shall limit ourself to the CP conserving case taking iPi = h, H 2 = H, and = A. 
And the couplings of the neutral Higgs bosons to the charged Higgs and W are: Sh = 
cos(/9 — a), Sh = — sin(/3 — a), and Pa = 1- Neglecting the contribution from the lightest 
Higgs boson h as in the decoupling limit cos(/3 — a) —)■ 0, we observe that in the high-energy 
limit the cross section of the subprocess behaves like 

aiwn ^ H+b) oc l 2 ci + Sngfi + PaQa]" + \SH9f1 - PaQaI" • (21) 

We note that the cross section suffers a huge cancellation between the top- and Higgs- 
mediated diagrams and a further cancellation between the Higgs-mediated diagrams. Taking 
the type H model as an example, we hnd 

a{W~^b -A- H~^b)\^ oc 4 tan ^(3 , 

Awn ^ = c(wn ^ oc 2 tan^ , 

AW^h ^ + ( 22 ) 

with Cl = —ShQh = —PaQa — tan/3. Note, for the W^b -A- H^b process, 

Awn ^ « o (V + 

while the high-energy behavior of the t-only, (t -|- i7)-only, and {t + A)-only amplitudes 
remains the same. 

Furthermore, independent of the type of 2HDMs we note that for W^b —?■ Pl^b (re¬ 
spectively W^b —)■ H^b) the s-channel (respectively the n-channel) top-exchange diagram 
interferes destructively with the t-channel Higgs-exchange W—A and W—H fusion diagrams. 
For demonstration we show in Fig. |^the cross sections for the subprocess W~^b -A- H^b (left) 
and W~^b — ?■ H~^b (right) as a function of center of mass energy ^/s in the MSSM. We illus¬ 
trate separately the top diagram alone, the sum of the top and pseudoscalar Higgs exchange 
diagrams, as well as all four diagrams. Note that “top-|-A” and “top-|-iJ” are extremely 
close to each other. It is clear from the plot that the dominant contribution is coming from 

^ Though we are working in the framework of 2HDMs, M^+, Mh, and Ma are very close to one another 
in the MSSM that will suit our purpose here. 
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vs (TeV) 


FIG. 3. The cross section as a function of center of mass energy for the subprocess W~^h —)• H^b 
(left) and W~^b —?■ H^b (right) in the MSSM for tan/3 = 30 and Ma = 400 GeV. 


the top diagram. It is also visible from the plot that the interference between s-channel top 
diagram and W-A fusion is destructive. The top contribution is reduced by a factor of 2 
by the W-A fusion diagram and same destructive interference takes place with W-H fusion 
diagram. We only show the sum of the top diagram and W-A fusion diagram in the hgure, 


that of the top and W-H fusion diagrams is almost the same. As expected from Eq. (22), 
in the case of W~^b —)■ W+6, after inclusion of all diagrams the total cross section drops by 
more than 3 orders of magnitude at large ^/s, as shown on the left panel of Fig. This in 
fact is due to the strong destructive interference of top diagram with the W-A and W-H 
fusion diagrams. Similarly, on the right panel in Fig. we illustrate the cross section for 
W^b —)■ H^b as a function of y/s. Again, as expected we can see destructive interference 
between w-channel top diagram and t-channel W-A and W-H fusion diagrams. We stress 
that the destructive interference in the fe-initiated process is less severe than the 6-initiated 
one, such that the total cross section for W~^b H^b is about one order of magnitude larger 
than that for W^b —)■ H^b. This is because the cancellation between the u- and f-channel 
diagrams is not as effective as in the s- and t-channel diagrams. For hxed and relatively 
small values of a/s, there are non-negligible parts of phase space in which the high-energy 
limits u/{u — ml) ~ 1 and t/{t — Mfj.) ^ 1 can not be achieved simultaneously due to 
the relation t + u = —s + Therefore, the 6-initiated process has an order of 

magnitude larger cross section. 
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FIG. 4. The pp —?■ H^bj cross sections in the 2HDM type II as functions of tan/3 (left panels) 
and charged Higgs mass (right panels) at LHC-14. The upper panels are the b initiated process, 
the middle panels are b initiated process while the lower panel are the sum of b and b. The charged 
conjugate panels are included in the plots. All panels are for the decoupling limit sin a = — cos/3. 
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B. For full process pp HHj 


In the previous subsection we have shown analytically and illustrated numerically the 
cancellation in the subprocesses W^h — hH^ and W^h —)■ H^h between the top diagram 
and W-A and W-H fusion diagrams using the effective W approximation. In Fig. we 
show the cross sections for the full 2 —)■ 3 processes qb —)■ q'H^b (upper panels), qb —)■ q'H^b 
(middle panels), and their sum (lower panels) as functions of tan (left panels) and Mh+ = 
Ma = Mh (right panels), including the charged-conjugate channels and after folding with 
the parton distribution functions Again, we separately show the contributions from the top 
diagram only, the top plus W-A fusion diagrams, the top plus W-H fusion diagrams, and all 
diagrams. We have assumed that we are in the decoupling limit sin a = — cos {3 and taking 
a spectrum of degenerate Higgs bosons Mh± = Mh = Ma as in the MSSM and the lightest 
CP-even Higgs boson h is the observed one with mh = 125.09 GeV. Thus, the diagram with 
h proportional to cos(/3 — a) does not contribute while the amplitudes associated with the A 
and H diagrams are the same up to a factor of 7 ^ in the Abb and Hbh vertex. Also, we can 
see that in the 6 -initiated subprocess (upper panels) the “top-|-A” curve completely overlaps 
with “top-l-if” curves but not exactly in the 6 -initiated one (middle panels). 

The upper panels in Fig. [^illustrate a very strong cancellation between the top diagram, 
and W-A and W-H fusion diagrams. Note the charged-conjugate channels qb —)■ q'H~b are 
included in it. On the other hand, the middle panels show a less severe cancellation between 
the top diagram, and W-A and W-H fusion diagrams where the charged-conjugate channels 
qb —)■ q'H~b are also included in it. Therefore, we still see that a strong cancellation occurs 
for the full process pp —)■ H^bj at the LHC-14, shown in the lower panels. Also, note that 
the total cross section is dominated by the 6 -initiated process when Mh+ > 200 GeV, where 
the internal top cannot be produced on-shell. 

It is clear from the left panels that the cross sections are enhanced for both small tan (3^1 
and large tan/?, the latter of which is associated with enhanced bottom Yukawa couplings. 
A dip indeed occurs around tan/3 6 in the tan/3 plots, which corresponds to where the 
top and bottom Yukawa couplings become similar m^tan /3 ^ mt /tan /3. We stress that our 
results are in good agreement with Ref. pj. 

^ Our numerical calculations of the several cross sections for the b- and 6-initiated full 2—^3 processes 
presented here are carried out by use of the Helicity Amplitude Method [23]. We compare our results for 

the total cross sections with those obtained using MadGraph [25] and find excellent agreements. 
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FIG. 5. The pp — )• H^bj cross sections at LHC-14 as a function of Ma in the 2HDM type I and 
III (left), and II and IV (right) for several values of tan/3. We have taken Mfj± = Mh = 300 GeV. 


As shown on the right panels, we emphasize that for Mh± ^ Trit — rrih the top quark can 
be produced on-shell as in single-top production and then decays into bH~^, Therefore, in 
the range of Mh± < mt — ruf, the top-exchange diagram completely dominates over other 
diagrams. On the other hand, for Mh± > rrit — mb the top quark is off-shell, and thus 
other diagrams also make signihcant contributions. Note that in the 6-initiated process 
qb —)■ q'H^b the top quark is never produced on-shell. 

It is well known that in the MSSM and for Ma > 200 GeV all the heavy Higgs bosons 
become degenerate Mh = Ma = Mh± and cos(/3 — a) —)■ 0. In general 2HDMs all Higgs 
boson masses are independent parameters. One can then identify the lightest CP-even with 
the observed 125 GeV Higgs boson and take the others Mh, Ma and Mh± as free parameters. 
In Fig. we show the total cross sections as a function of the CP-odd Higgs mass for a few 
values of tan /3 = 1 — 50 and Mh = Mh± = 300 GeV in 2HDMs types I and HI (left panel) 
and types H and IV (right panels). Note that for this choice of masses the production cross 
section is dominated by the 6-initiated process. For the values of the couplings in production 
cross sections, we refer to Table In types I and HI, the cross sections are insensitive to 
the CP-odd Higgs mass and they are suppressed by 1/tan^/3 with increasing tan/3. The 
largest value of the cross section is obtained at tan /3 = 1 and is of the order 27 fb. In types 
H and IV, one can see some sensitivity to the CP-odd Higgs mass. For Ma < 250 GeV, 
the cross section increases for lighter CP-odd Higgs mass and becomes almost constant for 
Ma > 250 GeV. The enhancement of the cross section for Ma < 250 GeV is amplihed with 
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FIG. 6. The transverse momentum pTf, distribution for pp —)• H^bj in 2HDM type II and IV 
at the LHC-14 for Ma = 100,250,400 GeV. We have taken the decoupling limit sin a = — cos/3, 
tan/3 = 30, and = Mh = 400 GeV. 

large values of tan 13. For Ma = 100 GeV and tan /3 = 30 one can reach a cross section of 
the order of 40 fb. 

In Fig. 1^ we plot the distribution of the b quark in the decoupling limit and for 
Mh± = Mh = 400 GeV and tan/3 = 30 for several values of Ma = 100,250 and 400 GeV. 
As one can see from the plot, the distribution is enhanced for light Ma = 100 GeV and 
for pti, < 200 GeV. The transverse momentum of the b quark is then a useful variable to 
separate the contributions between the top diagram and the W-A fusion diagram. One can 
require < 200 GeV to suppress the top-exchange contribution. Therefore, we can see 
that the W-A fusion diagram dominates for light Ma and at the lower pt^ region. 

C. Large tan/3 and LHC pp — )■ data 

At the LHG with 7 and 8 TeV, searches for the Higgs bosons <I> that decay into tan 
pairs, which in turn decay into those hnal states with one or two light leptons, have been 
performed [26l |27] for Higgs mass in the range [100, 900] GeV. Both ATLAS and GMS have 
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FIG. 7. Exclusion plots in the plane of {Ma, tan /3) for the 2HDMs type I, II, III and IV, using the 
ATLAS data for gg ^ A ^ . The excluded region is shown in red while the rest is allowed. 

some exclusion limits given as (j{gg —)■ $) x —)■ r+r“) as a function of the Higgs mass 

m$. These limits can be interpreted in the 2HDMs if we take the Higgs state $ as one 
of the neutral Higgs bosons of the 2HDMs: $ = h and/or H,A. In fact if h mimics the 
SM Higgs boson, a{gg h) x B{h ^ will not have any enhancement factor such 

as tan/3. Since we have observed that in order to enhance pp bH~^j cross sections one 
needs both non-degenerate A and H and also large tan/3, here we attempt to hnd what 
would be the largest possible value for tan /3 such that it is still consistent with data 

for 100 < Ma < 340 GeV and assuming that the heavy CP-even Higgs boson is rather 
heavy. A similar study with the 7 TeV data had been done in |2B] for 2HDM. Because of 
CP invariance the CP-odd Higgs boson A does not couple to WW or ZZ, and the partial 
decay widths into loop mediated gg, 77 , and 'jZ channels are highly suppressed. The decay 
channel A —)■ hZ, which is proportional to cos(/3 — a), will also be severely suppressed if 
we assume that /3 — a is close to the decoupling limit. Therefore, the CP-odd Higgs boson 
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predominantly decays into fermion pairs: qq, q = b,s, d, c, u and I = fi, e. 

In 2HDM-I and -IV, the coupling At^t~ is proportional to 1/tan/? while in 2IIDM-II 
and -III it is proportional to tan/?. On the other hand, from Table the coupling of Abb is 
proportional to tan/? in 2HDM-II and -IV but 1/tan/? in 2HDM-I and -III. Thus, it is clear 
that in 2IIDM-I (resp. II) both the production rate gg ^ A and the decay A —)■ are 

suppressed (resp. enhanced) for large tan/?. We then expect a strong exclusion for large 
tan /? in type II but not in type I. In all four 2IIDM types we expect some enhancement for 
small 0.5 < tan/? < 1 because Att is proportional to rrit/ tan/?. 

For a given CP-odd Higgs mass Ma and tan/?, the cross section of gg —)■ A, which 
only depends on these 2 parameters, is computed with help of SUSHI public code [29] . 
In the decoupling limit A —)■ Zh is vanishing, and if H —)■ ZH,tt} are closed, the 

branching ratio of H —>■ depends on tan /? and Ma only, and there is no sin a dependence. 

Therefore, the cross section gg ^ A ^ will depend only on tan/? and Ma- Hence, 

our exclusion from data can be given in the plane of (MA,tan/?). After computing 

the cross section gg ^ A times the branching fraction of A —)■ we compare our 

theoretical predictions with the ATLAS data [26]. Note that the ATLAS data were given 
for Ma = 90,100, ...340 GeV with steps of 10 GeV or even larger in some cases. Therefore, 
we have used linear interpolations for Ma values in-between the data. 

We draw in Fig. the exclusion region in {Ma, tan f3) plane for 2HDMs type I, H, HI, 
and IV, where Ma is in the range [100,340] GeV and 0.5 < tan/? < 50. In 2HDM types I 
and IV, there is no exclusion for tan/? > 1.5. The reason is that in type I, the production 
rate of gg ^ A and decay of A —)■ are both suppressed by 1/tan /?. Whereas in 2HDM 

type IV, the production rate is enhanced by the bottom Yukawa for large tan /? but the 
decay A —>■ r+r“ is suppressed by 1/tan/? which cancels the bottom enhancement in the 
production, which then gives no exclusion for large tan /?. In 2HDM type HI, similar to type 
I the production rate is suppressed by 1/tan/?. However, the branching fractions of A ^ qq 
are suppressed by 1/tan/? while B{A ^ MM) is enhanced for large tan/?. For this reason 
the exclusion in type HI is somewhat stronger than that in type 1. There is no exclusion 
for tan/? > 9.5 (resp. tan/? > 2) for Ma ~ 340 GeV (resp. for Ma ~ 112 GeV). On the 
other hand, the 2HDM-H receives both enhancement at large tan /? for the production rate 
gg ^ A and the B{A ^ Mt~). Thus, this gives a strong exclusion for tan/? > 22 for 
all Ma € [100,340] GeV. The low tan/? < 3 region (resp. tan/? < 1) is also excluded for 
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Ma = 340 GeV (resp. Ma = 150 GeV). 


IV. BEYOND TWO-HIGGS-DOUBLET MODELS 


Another interesting possibility to enhance the production of charged Higgs boson might 
be the case in which a 2HDM is not an ultraviolet (UV) complete theory and the UV cutoff 
locates far above the mass scale of heavy Higgs bosons In this case, taking one of the 
2HDMs as a low-energy reference model, the relevant interactions may be parameterized as 


^Hibb — 


^H±tb — d 


b gf + i gf 75 ) b Hi , 

b{^LCLPL + ^RCRPR)tH~ + h.c., 


2 mw 

grrib 


\/2mw 

= —^{^sSi + i^p Pi) 


H- 




+ h.c.. 


(24) 


In the MSSM, for example, including the tan /3-enhanced SUSY threshold corrections to the 
down-type Yukawa couplings, we have 

1 


ef = er = a = 


I + Kb tan (3 ’ 


^R = ^s = ^p = l, 


(26) 


with 


^6 — G + 

where eg and ep are the contributions from the sbottom-gluino exchange diagram and from 
stop-Higgsino diagram, respectively. Their explicit expressions are 


2a. 


G = l^sl ), 


Ch = 


\ht 






where M 3 is the gluino mass parameter, ht and are the top-quark Yukawa and trilinear 
couplings, respectively. 

Without loss of generality we choose the 2HDM type H as the reference model, and show 


S P 

the change in production cross sections with the variations in the couplings ^i ’ and 


We used Eq. (25) as the guidance. We hrst show the ratio of the cross sections for varying 
between —2 and -|-2 to the cross section at the 2HDM type H values, i.e., 
= if = ii = 1 in Fig. ^ for tan/3 = 1,30. The tan/3 = 1 curves show almost no 


We note that this kind of enhancement arising from unitarity violation might be dangerous and should be 
taken with caution, because a UV-complete model might contain new interactions to restore the unitarity. 
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FIG. 8. Ratios of cross sections for varying to the cross section at the 2HDM type 

II value = 1), showing the enhancement due to the modification of the 2HDM as 

suggested in the text. 
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FIG. 9. Ratios of cross section with varying (left panel) or (right panel) to the cross section 
at the 2HDM value (i.e. = 1 or = 1 respectively) by keeping all other parameters at their 

2HDM type II values. 


sensitivity to because the process is dominated by the top-Yukawa term. On 

the other hand, the tan/3 = 30 curves are dominated by the bottom-Yukawa term. It is 
obvious that the ratio is close to zero for if = if = = 0, and is one for if = if = ii = 

The ratio grows as the square of the couplings around 0 to almost 4 at if = if = ii = ±2. 

If for some higher scale dynamics such that if and if do not change in the same manner, 
we show the effects on the cross sections in Fig. On the left panel, we show the ratio of 
cross sections with varying if between —2 and +2 to the cross section at the 2HDM value 
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(i.e. = 1) by keeping all other parameters at their 2HDM type II values. Again, the 

sensitivity at tan/? = 1 is negligible, while it becomes quite nontrivial for large tan/? = 30. 
As we have shown in Sec. Ill that the W-A diagram interferes destructively with the top 
diagram, we can now turn the destructive interference into constructive one by reversing 


the sign of Furthermore, when is negative the second term in Eq. (21) would not 
vanish. It is then very clear that the ratio becomes quite large at negative At = 0, 
the ratio is already larger than 1 because no interference comes from the W-A diagram. 
Similar behavior occurs for as shown on the right panel in Fig. The ratios that can 
attain are very similar to those of 


V. CONCLUSIONS 


We have performed the study of b- and 6-initiated processes of pp —)■ jH^b/b in the 
2HDM framework at the LHC-14 in the decoupling limit (sin a = — cos /?), which is favored 
by the current Higgs data. We have identihed strong cancellations between the top diagram 
and the W-Hi{Hi = h,H,A) diagrams which rendered the process very suppressed. The 
cancellation is indeed the strongest when the A and H are degenerate and in the decoupling 
limit. We pointed out that if the pseudoscalar Higgs boson A is much lighter than the 
CP-even Higgs boson H, the cross section of charged-Higgs production can be substantially 
enhanced, because the cancellation is no longer complete. We have explicitly obtained the 
exclusion in parameter space of {Ma, tan /?) for 2HDM types I to IV based on the LHC data 
on a{g —)■ <F) x i?(<F —)■ r+r“). In the allowed paramete space, the size of production cross 
section can be as large as 0(50) fb for Ma = 100 GeV and tan/? = 30 for types H and IV. 
This is the main result of the work. 

We offer the following comments on the Endings of this work as follows. 


1 . The 6-initiated process for production of in qb —)■ qH^b suEers from a very 
strong cancellation between the top diagram and the W-Hi diagrams. However, the b- 
initiated process for production of in qb —)■ qHM suEers a less severe cancellation, 
mainly because of the w-channel top-exchange instead of s-channel. 


2. The strong cancellation is dictated by the absence of the unitarity-breaking terms and 


we find the sum rules expressed by the relevant Higgs couplings, see Eq. (17). 
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3. For Mh± < rrit — rrib the top diagram completely dominates as the top quark is 
produced on-shell as like that in single-top production. However, when Mh± > nit—nib 
the W-Hi fusion diagrams also contribute. 

4. In the future study, we shall make use of the special kinematics, e.g, the pt^ distri¬ 
bution, to discriminate between the top and the W-Hi fusion diagrams. The goal is 
to isolate the effect of light pseudoscalar Higgs boson, which is still allowed by the 
current data. 

5. Current LHC data on a{gg —)■ <h) x H(<h —)■ constrains the parameters of the 

2HDMs. We found that the data constrains the most on type H because both the 
production process and the decay are tan/? enhanced. Yet, there are sizable allowed 
parameter space between tan/? = 3 and 22. The second most constrained is type HI 
because the production rate is suppressed by 1/ tan (3 but the decay is enhanced by 
tan 13. Types I and IV have the most available parameter space. 

6 . The process pp —)■ jH^b/b that we consider in this work would be more interesting 

for type H and IV because of larger cross sections. Especially type IV has the least 
restriction from the current LHC data on a{gg —)■ <h) x H(<F —)■ and it can 

allow cross sections as large as 0(100 — 300) fb for tan /3 = 30 — 50 and Ma = 50 — 100 
GeV. 

7. When the 2HDM is a low-energy limit of some ultraviolet (UV) models, the integrity 
of the Yukawa couplings may change. For example, in the MSSM the SUSY particles 
can largely change the bottom-Yukawa couplings with strong and weak interactions. 
Varying the bottom Yukawa coupling of either A or H gives non-trivial behavior for 
the production cross sections. 
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